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Abstract 

Recent numerical work on the gravitational collapse of a 5D (4+1) Yang-Mills instan- 
ton has provided numerical evidence that the free energy F = E — TS = E/3 of a 5D 
Schwarzschild black hole of mass E can be obtained classically via the Lagrangian. Al- 
though there is no Hawking radiation, these numerical results suggest that the quantity 
TS has a classical meaning. We investigate this association for the physically relevant case 
of 3+1 dimensional collapse. We track numerically the negative of the total Lagrangian -L 
during the gravitational collapse of a massless scalar field to a Schwarzschild black hole in 
isotropic coordinates. Wc show that — L approaches the free energy F = E — TS = E/2 of 
a 4D Schwarzschild black hole to within 5%. We also show that the matter contribution 
to the free energy tends towards zero so that the entropy at late stages of the collapse is 
gravitational in origin. The entropy S makes a negative contribution to the free energy 
and this feature is observed in our numerical simulation. There is a pronounced dip (neg- 
ative contribution) in a thin slice just inside the event horizon precisely where the metric 
field is nonstationary. This is in accord with recent work suggesting black hole entropy is 
connected with the nonstationary phase space hidden behind the event horizon. We also 
obtain thermodynamic results for the 5D collapse of a massless scalar field which confirms 
that previous 5D results are universal and independent of the type of matter undergoing 
the collapse. 
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1 Introduction 



The study of thermodynamics during gravitational collapse provides a window into the process 
by which black holes reach their final thermal state. The majority of the work to date has 
consisted of analytical and semi-analytical studies of black hole entropy and Hawking radia- 
tion during shell or dust collapse [TJ- [8] , and the approaches have been quantum mechanical 
or semi-classical. Recently, a classical approach, based on the free energy, was carried out [9]. 
The authors conducted a numerical study of classical thermodynamics during the gravitational 
collapse of a 5D (4 + 1) Yang Mills instanton to a Schwarzschild black hole in isotropic coor- 
dinate^]. The function F = —dl/dt = —L was tracked during the collapse process where / is 
the total gravitational plus matter action, L is the total Lagrangian and t is the time mea- 
sured by a stationary clock at infinity. At late stages of the collapse, F approaches a constant 
that can be identified with the free energy F = E — TS where E is the ADM mass of the 
Schwarzschild black hole and S and T are the entropy and Hawking temperature respectively 
(the identification with the free energy is discussed later). The entropy and temperature for 
a 5D Schwarzschild black hole are known from standard black hole thermodynamics [114 [12] : 
T = %/{Att ro), S = 4tt 2 ^/(HGs) where tq is the gravitational radius given by r 2 = 2G$ E/(3tt) 
and G5 is Newton's constant in 5D. The product ST is equal to 2E/3 so that the free energy 
in 5.D is equal to E/3. The authors showed that at late stages of the collapse the function 
F = —L approached a numerical value close to E/3. Although there is no Hawking radiation, 
the numerical results suggest that the product TS can be obtained classically via the gravita- 
tional Lagrangian. Note that though S and T each contain h, the product S T appearing in 
the free energy does not. 

We investigate the association between — L and the free energy E — TS for the physically 
relevant case of 3+1 dimensional collapse (for completeness, we also consider a massless scalar 
field in 4+1 dimensions). We track numerically the function — L during the classical gravita- 
tional collapse of a AD and 5D massless scalar field to a Schwarzschild black hole in isotropic 
coordinates. As far as we know, this is the first time that the negative of the Lagrangian has 
been tracked numerically during 3+1 dimensional collapse. The entropy and temperature of 
the 4D Schwarzschild black hole of mass E is well known and given by S = 4irG E 2 /{h) and 
T = h/(8irGE) PUE2]. The free energy F = E—T S is therefore equal to E/2. At late stages of 
the collapse, when thermal equilibrium has been reached, —L is observed to reach a constant 
numerical value that is within 5% of the free energy E/2. We also study the gravitational 
collapse of a 5D massless scalar field. This case differs from the 5D Yang-Mills instanton of 
[9] because a static massless scalar field may undergo gravitational collapse to a black hole 
or disperse depending on its initial configuration [13] . In our simulation, this is governed by 
the value of the scale parameter A appearing in the initial state for the scalar field. For the 

1 The 5D Yang Mills instanton acts like dust [10] so that any initial static state is guaranteed to undergo 
gravitational collapse to a black hole. 
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case of collapse to a black hole, we obtain a value for —L which is again within 5% of the free 
energy (E/3 in 5D). This confirms that the previous 5D thermodynamic results obtained for 
the collapse of a Yang-Mills instanton [9] are independent of the type of matter undergoing the 
collapse. Our results together with previous work [9] provides considerable numerical evidence 
for associating a purely classical quantity, the negative of the gravitational Lagrangian, to the 
free energy E — TS of a black hole. 

We show analytically that the exterior region of a Schwarzschild black hole by itself makes a 
contribution E to —L. It follows that the interior makes a negative contribution (from -TS). 
This feature is observed in both our AD and 5D simulations. In a thin region just behind the 
event horizon there is a significant dip (negative contribution) to the function —L. Moreover, 
one observes that the metric fields are nonstationary in this region (and are static everywhere 
else). This is in agreement with recent work [14] that suggests that black hole entropy is 
connected with the nonstationary phase space hidden behind the event horizon. It is also in 
agreement with the entropy obtained via explicit counting of internal microstates within the 
context of canonical quantization of the Schwarzschild [K] and the Reissner-Nordstrom black 
hole [16] . That the entropy of a black hole is connected to its interior is in line with Bekenstein's 
original view that black hole entropy is a measure of an outside observer's ignorance of the 
internal configurations hidden behind the event horizon [11] , 

1.1 Metric in isotropic coordinates 

We work in isotropic coordinates, where a spherically symmetric time-dependent 4L> metric 
takes the form 

ds 2 = -N(r,t) 2 dt 2 + ^{r,tf{dr 2 + r 2 dfl 2 ) . (1) 

TV is called the lapse function and we refer to ip as the conformal factor. The coordinate time 
t is the time measured by a stationary clock at infinity. During our simulation of gravitational 
collapse, N and ip remain finite and continuous everywhere and the metric (pQ) covers the entire 
spacetime (inside and outside the event horizon). We assume asymptotic flatness so that N 
and ip approach unity at infinity. Black hole formation will be interpreted to coincide with the 
formation of an apparent horizon. This occurs when N crosses zero; light cannot reach a point 
where iV = from any other point in a finite time 117] . 

The metric in the exterior region of the Schwarzschild black hole in isotropic coordinates is 

pairs] 

l2 (l-GM/2r) 2 l2 / GM\4 

The event horizon is located at r = GM/2. The region r > GM/2 of metric ([2]) covers the 
exterior region of the Schwarzschild black hole. However, the region r < GM/2 covers the 
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exterior region a second time, not the interior, so that metric ([2]) represents a double covering 
of the static exterior region. Metric ([2]) is obtained via the coordinate transformation r' = 
r (1 + GM/2r) 2 where r' is the radius that appears in the standard Schwarzschild form ds 2 = 
-A(r') dt 2 + B(r') dr' 2 + r' 2 dVL 2 with A(r') = 1 - 2GM/r' and B(r')= A(r'y l . The minimum 
value of r' is 2GM so that the interior region r' < 2GM is not covered. This is to be expected 
because the metric ([2]) has a timelike hypersurface-orthogonal Killing vector in the region 
r < GM/2 and is therefore static in that region. In contrast, all the Killing vectors in the 
interior region of the Schwarzschild black hole are spacelike and the region is nonstationary. 
To cover the interior region of the Schwarzschild black hole in isotropic coordinates, the metric 
coefficients in the region r<GM/2 need to be time-dependent |18] , 

At late stages of the collapse process, the metric ([pj approaches the static form (J5]) only in 
the exterior region r> GM/2. In the interior region (r <GM/2), the metric functions N(r,t) 
and 4>{ r it) are time-dependent at late times in accord with the nonstationary nature of the 
Schwarzschild interior. The areal radius is R = ip 2 r and the Schwarzschild central singularity at 
R = corresponds in isotropic coordinates to ip{r,t) — > in the interior region as t — >■ oo. The 
function tjj(r, t) is continuous everywhere so that the interior and exterior values must match 
at the event horizon r = GM/2 (where the analytical solution yields ip = 2). Therefore, as one 
crosses the event horizon from the exterior towards the interior, ip must decrease from a value 
of 2 at the event horizon to a value approaching zero inside. This is observed in our numerical 
simulation. At late times, ip has a small constant value throughout the interior except for a 
small region just inside the event horizon where it increases sharply to a value close to 2. In 
this small region, ip is nonstationary whereas it is basically static everywhere else. We will 
see that the entropy of the black hole originates from the nonstationary region just inside the 
event horizon. 

1.2 The free energy 

The thermodynamic function tracked during our simulation is F(t) = —dl/dt = —L. It was 
argued in [9] that if F(t) approaches a constant at late times in the collapse process, then 
it can be identified with the free energy F = E — TS of the black hole. The argument is as 
follows: if F(t) approaches a constant at late times, then the action I becomes linear in time 
and I = —F (t — to) in the time interval between to and t. The corresponding Euclidean action 
Ie is obtained from the Lorentzian action / via the relation lE = Ht->ir = F (r — tq) where r is 
the Euclidean time. The quantum partition function Z is equal to the Euclidean path integral 
[20| with time interval r — tq = 1/T where T is the temperature and with fields (j> obeying 
the periodic boundary condition 4>{tq) = 0(ro + 1/T). For a single (classical) trajectory, the 
partition function reduces to Z = e~ lE = e~ F l T . We see that F is indeed the free energy 
F = —TlnZ. The negative of the Lagrangian is therefore identified with the free energy 
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F = E — TS at late times. This identification was used in the numerical work of [9] where 
—L at late times was compared to F = E/3 for the case of spherical collapse of a Yang Mills 
instanton to a 5D Schwarzschild black hole of mass E. 



2 4D thermodynamics during classical gravitational collapse 
2.1 Equations of motion in isotropic coordinates 

In this section we obtain the equations of motion for the metric and matter fields in isotropic 
coordinates. We consider a massless scalar field x coupled to gravity in AD. The Klein-Gordon 
action for a massless scalar field is given by 

Ikg = -\ j V^QST^XdyX d A x = J ^C KG d 4 x (3) 

where Ckg = (—^/2)g fJ ' u d^x9uX- The action for gravity is Iq = I eh + Ib where I eh = 
16 * G f ^J—g R d 4 x is the Einstein-Hilbert action and Ib is a boundary term that is required if 
the variation of the total action I = Iq + Ikg with respect to the metric is to reproduce the 
Einstein field equations |21j . We are assuming spherical symmetry so that x = x( r >^)- 

Variation of the action / with respect to the metric leads to the Einstein equations, G^ v = 
8tt G T^ u where G^ v is the Einstein tensor and T^ u is the energy-momentum tensor given by 



(dC K G\ 



)+C K G9„u- (4) 



The non-zero components of and G^ v for the isotropic metric ([I]) are 



, N± I XI , 2^ ) T - r f 2C + K I T , — T, — yy' 

1 tt — o \ t TTI I , J-rr — o Jjl ) , J-rt — J-tr — XX 



12 , x 2 \ . T _ i^(x 12 

r« = ; T H = sinHT 66 

G tt = ^(3r?/>V - 2^'N 2 - r^'N 2 ) 

G rr = v ^ w (2r^' 2 N 3 + 2ip'N 3 ?p - 4np 2 N^ + 2rN^ 5 - 2r^Nip 5 + 2rN'ip'N 2 ip + N'N 2 ip 2 ) 

G rt = G tr = ^{^'N -ijj'NiP + ipN'iP) ; G H = S m 2 9G e e 

Gee = ^1 (-2^'N^ - 2r^"N 3 ^ + 8r^ 2 ^ 4 

+2rij' 2 N 3 - irNipip 5 + ^rNtjj 5 ^ - N'N 2 ifj 2 - rN"N 2 ip 2 ) . 
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2.1.1 Energy and momentum constraint and evolution equation for tjj 

The Einstein equation Gtt = k 2 Tu where k 2 = 8ttG yields 

-L(3^V - - r^'N 2 ) = ^fl$ + ^). (5) 

The above equation is the energy constraint and it can be expressed in the convenient form 

-±V^ = A-f. (0) 



In (|6]), V 2 ^ is the flat space Laplacian given by 

VV=^V 2 ) = ^ + /, (7) 
and K is the trace of the extrinsic curvature given by [21] 

K = h*K* = ^h*d t h ab = -^ (8) 

where h ab is the three-metric induced on a spacelike hypersurface at an instant of time t. 
The term £ is interpreted as the energy density of the matter field \'- 



The evolution equation for the field ip is given by ([8] 

N ~ 6~~ 



(10) 



If K, N and tp are known at time t, ip can be evaluated at the next time step. The energy 
constraint equation ([6]) is not an evolution equation. It is used to monitor the accuracy of the 
simulation. 

The Einstein equation G r t = k 2 T T t yields the momentum constraint equation 

4 



ip 2 N 

Using (|10p . this can be expressed in the more compact form 



(ip^'N -ip'Nip + ipN'tp) = k 2 xx'- (11) 



- = -^x'- (12) 
3 2 N X 1 ' 
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2.1.2 The evolution equation for K 



The Einstein equation G rr — k T rr yields 
2 



( 2 r V' 2 iV 3 - 4r-0 2 + 2?//iV 3 V> + 2rNip^ 5 - 2r^iV<0 5 + 2rN'ij)'N V + N'N 2 ^ 2 



rtfj 2 N 3 



2 V^+^J' < 13 » 

Using the expression for the energy density £ given by ([9]) and the expression for K given by 
(|10p . the above equation reduces to a time evolution equation for K, 

K K 2 3 2 „ ^j'/V 7 ' 1\ JV' 1\ , x 

iV = T" + 2 KS + (T + r ) • (M) 



Following a prescription by Finelli and Khlebnikov [T7] , one can improve the numerical stability 
of equation (fT4"j) by adding to it the term — |w(k 2 £ — ^ +4^^)) where w; is a constant. The 
added term is equal to zero because of © and therefore does not alter the equation. The final 
evolution equation for K is 

iV = — (1 + " } + 2 K ^ (1 " W) ~ 6 T 5 ^ + r)- 3 N4A \— + - r )- 6W ~4F ■ (15) 
In our simulation, w is set to unity. 



2.1.3 Lapse function N: ordinary differential equation 

There are two metric functions that appear in the isotropic metric dU): ip(r,t) and N(r,t). 
The function ip{r, t) is a dynamical variable and its evolution is determined by (jlOp with K 
evolving via (|15p . The lapse function N(r, t) is not a dynamical variable and there is no explicit 
evolution equation for it. Nonetheless, each time step yields a new spacelike hypersurface St 
and one must readjust the value of N(r,i) at each time step. This is governed by the Einstein 
equation Gqq = k 2 Tqq i.e. 

= 1 = ( -2ib'N 3 ib - 2r^"iV 3 ^ + 8rV> 2 A^ 4 + 2r^ /2 iV 3 
ip l N 6 r 

.2/4 /2 -2 

-ArNtptp 5 + ArNtp 5 ^ - N' N 2 ^ 2 - rN"N 2 ip 2 ) = - . (16) 

2 V ip^ A'' 2 / 

The sum of the terms with time-derivatives can be replaced by spatial derivatives via equation 
(|13p . The lapse function N is then governed by an ordinary differential equation containing 
only derivatives with respect to the variable r, 

V^M^) + N dr (v^) = • (17) 
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2.1.4 Evolution equations for the matter field 



Variation of the action Ikg given by ([3]) with respect to the massless scalar field x yields the 
Klein-Gordon equation in curved spacetime, 

9^^^uX = 0^d^(^g^d l/X ) = 0. (18) 
For the metric ([1]), ^J—g = N ij)® r 2 sin((9) and the above equation becomes 

±d r (Nip 2 r 2 X ') = d t ^ 6 x/N). (19) 
Defining a new function, p(r, t) as 

P = V> 6 ! (20) 



equation (|T9[) can be expressed as 

The evolution equation for X is given by (|20p with the function p(r, t) evolving according to 
(|2ip . We now express the energy density ([9]) and the momentum constraint equation (112h in 
terms of p, i.e. 

2 V t/; 4 ^ 

and 



K' K 2 i 

Y = ^ xp - (23) 

From the form of (|22p . we see that p acts like a momentum conjugate to the field x- 



To summarize, the evolution equations for tp and K governing the metric field are given by 
()10p and (I15p respectively whereas the evolution equations for x an d p governing the matter 
field are given by (|20p and (|2ip respectively. The lapse function N does not obey an evolution 
equation but is evaluated at every time step via the ordinary differential equation (|17l) . Once 
the boundary conditions and initial states are specified, the time evolution of all the fields is 
unique. 

2.2 Expression for the ADM mass 

During the gravitational collapse process the total energy (the ADM mass) remains constant 
and this allows us to monitor the accuracy of the simulation. We therefore need to derive an 
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expression for the ADM mass in isotropic coordinates (pQ) . The ADM mass of an asymptotically 
flat spacetime is defined as |2T] 



Madm = —t- hm f {k- k )^d 2 0, (24) 



where St is the two-sphere at spatial infinity, oab is the metric on St, k = a AB kAB is the trace 
of the extrinsic curvature of St embedded in Et, the three-dimensional spacelike hypersurface 
at constant t, and ko is the trace of the extrinsic curvature of St embedded in flat spacetime. 
The trace k is given by k = V a r a [21] where r a is the unit vector normal to the two-sphere 
boundary St- In AD isotropic coordinates, a quick calculation yields k = Aif)' / ^ + 2/(ip 2 r), 
k Q = 2/(V> 2 r) so that k - k =4^'/^ 3 . With [fad 2 Q = iji A r 2 sin 6 d6 <fy, the ADM mass 
reduces to 

2r 2 2R 2 
Madm = —gTp'ip\ r=R = —g-ip'\ r=R (25) 

where R is assumed large (infinite limit) and we used the fact that ip — >■ 1 as R — > oo. Using 
the energy constraint equation ([6]), one can express the ADM mass ()25[) in the integral form 

Madm = E tot = ^ J* (s - ^) ^ r 2 dr (26) 

where K and <f are given by ([8]) and ([9]) respectively. The ADM mass represents the total 
energy Etot of matter plus gravitation. During the collapse process, it remains a constant. The 
black hole mass E is obtained via the relation ro = GE/2, where ro is the location where the 
lapse function N crosses zero at late times. The black hole mass E will be less than the ADM 
mass Etot because in the numerical simulation part of the total energy appears in an outgoing 
matter wave. 



2.3 Initial states and boundary conditions 

We choose a static initial state, where K and p are zero. From (jlOp and (|20p . this implies that 
ip = x = at time t = 0. For this initial state, the momentum constraint (|23p is automatically 
satisfied. The initial field configuration for the scalar field \ is chosen to be 

8A 2 r 4 

X(t = 0,r)= (A2 + r2)4 , (27) 

where A is a scale parameter, x decreases rapidly to zero at the two extreme limits (i.e. as 
r — > and r — > oo), reaches a maximum at r = A and has a half-width of ~ 1.3 A. The 
initial energy density is concentrated in a shell of approximately this width and its peak value 
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increases as A decreases. For A sufficiently small, the self-gravitational pull becomes large 
enough to initiate a gravitational collapse to a black hole. 

The initial state for the conformal factor ip is obtained by solving the energy constraint ([6]). 
Asymptotically, the spacetime is flat: •0—7-1 and N — > 1. To solve for ip with this boundary 
condition, we express it as a power expansion in the parameter Q = k 2 i.e. 

■ip = 1 + C0i + C 2 ^2 + • • • • (28) 

The actual expansion parameter turns out to be £ divided by some power of the scale parameter 
A. By substituting the energy density £ given by (|2"2"j) and the \D flat-space Laplacian ([7]), 
along with the static assumption (K = p = 0) into the energy constraint equation ©, one 
obtains 

^r(rV) = 4^' ( 29 ) 
r z 8 

Substituting the expansion (|28|) into (|29p and matching the left hand side with the right hand 
side order by order in £ we obtain 



+ 



A 3 40320 r (r 2 + A 2 ) 8 
C 2 1 



1575 r 15 A + 12075 r 13 A 3 + . . . + 1575 (r 2 + A 2 ) 8 tan _1 (r/A) 



2010133125tt 2 - — ^ — r^l - r A 2 (2967339375 r 30 

r(r z + A^) 10 



A 8 10538886758400 
+ 47828405625 r 28 A 2 + . . . + 60640845881 r 2 A 28 + 4113561991 A 30 ) 
+ 1276275 A(r 2 + A 2 ) 8 (825r 16 + 4500 r 14 A 2 + . . . + 15436 r 2 A 14 + 2717A 16 ) tan^r/A) 

+ 4020266250r (r 2 + A 2 ) 16 [tan~ 1 (r/A)] 2N 



+ ... 

While we only show here the first two order terms in the expansion, our results are obtained to 
six orders. The initial state for N is obtained numerically by solving its associated ODE (|17p 
for the initial values of ip and x defined above. Starting at the outer boundary r = R, where 
the spacetime is flat and N = l, we iterate backwards to obtain N. 

Boundary conditions need to be imposed at r = to ensure regularity of the solution. To ensure 
asymptotic flatness, we also impose boundary conditions at "infinity" , the computational outer 
boundary r = R. The regularity conditions are x'(0,t) = and K'(0,t) = 0. The boundary 
conditions imposed at the outer boundary are N(R,t) = 1, K'(R,t) = and p'(R,t) = 0. 
Together with the initial matter and metric states, these lead to a unique evolution. 
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2.4 Expression for the total Lagrangian in isotropic coordinates 



We track the function F = —L where the total Lagrangian L is a sum of the Klein-Gordon and 
gravitational Lagrangian. The Klein-Gordon (KG) matter Lagrangian is given by 



Lkg = j V^g£-KGd 3 x = ~ j V^gg^ w d^xduXd 3 x 

r R /v /2 v 2 \ ( 30 ) 

Jo V W 4 N*J 



where R is the radius at the outer boundary. The gravitational action Iq is a sum of the 
Einstein-Hilbert action plus a boundary term Iq and a nondynamical term Iq [21| : 

I G = —— j ,/^jRd 3 x dt + I B + I = j Lq dt (31) 



16vrG . 

where Lq is the gravitational Lagrangian and contains only first derivatives of the metric 
functions. The Ricci scalar in isotropic coordinates is given by 

R = ^ (l8n/> 2 iV^ 3 + 6r^V 4 -8^'N 3 -4rip"N 3 (32) 

rN A ip b 

- GrNipip 4 - rN"N 2 ip - 2rN'ij/N 2 - 2N'N 2 ip) . 

Lq is obtained by integrating out by parts the second derivative terms appearing in R [9]. 
This yields 

L G = — %r 2 i>'N'i> + 8r 2 NiP 12 dr. (33) 

4G Jo \ TV / 

The total Lagrangian, L = Lq + Lkg, is finally given by the expression 

L = _L f R ( Rr 2 ,// N' ih + Rr 2 N ,// 2 - 24 ^ \dr - 9. tt N ^ r 2 (?£ - 2^ 



{sr 2 ^N'^ + Sr 2 N^' 2 - ^ ^ * ) dr -2tt ^ V - |-) . (34) 



2.4.1 Static exterior and nonstationary interior contributions to the free energy 

We now show that for a static scenario where ip = x = everywhere, -F = — L reduces to the 
mass E so that the product T S appearing in F = E — T S is zero. Setting ip and x to zero in 
the total Lagrangian (f34|) yields 

-Static = - ^ J* (8r 2 xj}' N' i/> + 8 r 2 N 2 ) dr + 2 vr jf * TV V> 2 r 2 X ' 2 

= -%T$ N i\n+Tri I (l6rip'Nip + 8r 2 ip" Nip)dr + 2ir [ N ip 2 r 2 x' 2 dr (35) 
^ 4G Jo ^ 'jo 

2R 2 „ 



G 
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where we integrated by parts and used the energy constraint ([6|) with ift = x = 0. The last 
line is simply the expression for the ADM mass (note that asymptotic flatness implies that 
N = tl) = l at r = R). 



The Schwarzschild black hole is not static everywhere: it is not static in the interior region, 
inside the event horizon. Its free energy will therefore not reduce to E as above. This agrees 
with the fact that the Schwarzschild black hole has a nonzero temperature T and a nonzero 
entropy S so that T S ^ 0. Nonetheless, we can still show that the static exterior region 
(r > GM/2) contributes by itself E to the free energy. It then follows that the — T S contribution 
stems from the nonstationary interior. The metric (J2]), which describes the exterior region of 
the Schwarzschild black hole, corresponds to a conformal factor of ^(r) = 1 + GM/2r and lapse 
function N(r) = (1 — GM/2r) /(l + GM/2r). The Schwarzschild exterior is a vacuum where x' 
and x are zero. Substituting these values into the Lagrangian (|34p and integrating from the 
event horizon at r = GM/2 to infinity yields 



where M and E are the same quantities (we work in units where c = 1). The free energy 
F = E — T S of a Schwarzschild black hole separates nicely into two contributions: a positive 
contribution E from the static exterior region and a negative contribution —TS from the 
nonstationary interioJl. This feature of a Schwarzschild black hole is observed in our numerical 
simulation. At late stages of the collapse, the metric field tjj is nonstationary only in a thin slice 
just behind the event horizon. This small nonstationary region makes a negative contribution 
to the free energy (a large dip is observed just behind the event horizon in the accumulation 
plot of the free energy). The exterior is basically static and makes a positive contribution. 

That black hole entropy stems from the nonstationary interior region is discussed in ref. |14j . 
Black hole entropy is a measure of an outside observer's ignorance of the internal configurations 
hidden behind the event horizon [11]. In [T3], the internal configurations were identified as 
points in phase space, the classical microstates [h a f,,P ab ] where h a b is the three-metric in the 
usual 3 + 1 decomposition [21j and P ab is its momentum conjugate. A nonstationary interior 
implies that h a b is time-dependent and that there is a continuous set of classical microstates 
[hab(t), P ab (t)] parametrized by the time t. The outside observer cannot track the parameter t 
and is clearly ignorant of which classical microstate the black hole interior is in. This implies 
a non-zero black hole entropy. 

2 It should in principle be possible to calculate analytically the contribution of the interior. However, this 
would be a lengthy calculation more appropriate for a separate analytical study of the free energy of stationary 
black holes. 




(36) 
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2.5 4D numerical results 



We work in geometrized units where G = c=l. In these units, k 2 = 8tt and mass, energy, time 
and distance have dimensions of length. In our simulation, these quantities are expressed in 
units of the ADM mass Etot which is normalized to unity. The scale parameter used in our 
simulation is A = 3/2 (the actual value is not that important as long as it is small enough to 
lead to black hole formation instead of dispersion which is the case here). It is convenient to 
make the change of variable r — > x = pr^. This reduces the overall number of grid points 
while maintaining a high density of points in the region where large gradients are present at 
late times. The numerical code is based on a fourth order Adams-Bashforth-Moulton explicit 
scheme |22[ 123] with space and time increments equal to Ax = A£ = l x 10 -4 . The total energy 
Etot is calculated from the integral form of the ADM mass (|26p . As already mentioned, the 
mass E of the black hole will be less than the total energy Etot because some energy appears 
in an outgoing matter wave. The ADM mass Etot is a conserved quantity which should remain 
constant during the evolution and we therefore use it to monitor the accuracy of the simulation. 
We show the plots up to t = 20 at which time the function F has plateaued and Etot is just 
beginning to deviate from unity. 

The metric functions ip(r,t) and N(r,t) as well as the matter function x( r ;0 are plotted as 
functions of r (on a log scale) at different times t in figures [U [2] and [3] respectively. At late 
times, in the exterior region, the conformal factor ifi (fig. [1]) can be seen to match closely the 
analytical form ip = 1 + M/2r of the Schwarzschild metric ([2]). This provides a check on the 
numerical simulation and confirms that the collapse process has indeed led to a Schwarzschild 
black hole. The event horizon is at r = 0.334, the location where the lapse function N(r, t) (fig. 
[2]) crosses zero at late times. After N crosses zero, it is negative throughout the interior region. 
Inside, it approaches zero from below but never actually crosses zero. The spacetime is well 
behaved throughout the evolution and all curvature scalars are finite. From the plot of the 
matter field x (fig- ElK on e observes that at t = 4, the initial shell configuration has expanded 
into the interior to form a spherical ball of matter. By that time, it has not yet collapsed to a 
black hole since N has not crossed zero yet ( N crosses zero for the first time at around t = 6) . 
At late times, the scalar field x nas collapsed to a thin shell [21] near the event horizon except 
for an outgoing matter wave in the exterior region. The crest of the wave at t = 8 and r ~ 5 
propagates outwards and reaches r ~ 10 at t = 14 and r ~ 15 at i = 20. 

In the function F = —L, there are two contributions to the total Lagrangian L given by (|34D : 
one from gravitation (Lq) and one from matter (Lkg)- Both Lq and Lkg ar e expressed as 
integrals over r where r ranges from to the outer computational boundary R. It is instructive 
to see where F accumulates as the integral ranges from to R. These accumulation plots are 
shown at different times t for gravitation (fig. H]) and matter (fig. [5]). At late stages of the 
collapse, the matter contribution to the free energy is almost zero: the free energy stems almost 
entirely from the gravitational part and we see that black hole entropy is gravitational entropy. 
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In fig. [31 at the late time t = 20, there is a dip in a thin region K of thickness Ar = e just 
inside the event horizon (0.334 — e<r< 0.334). The region 5? is nonstationary: as can be seen 
from fig. [6j tp peaks in Ji, is zero in the exterior region and tends with time towards smaller 
values in the interior region outside K. In short, tp is basically static everywhere except in the 
nonstationary region 5R. In the previous section we saw that the nonstationary interior makes 
a negative contribution to the free energy of a Schwarzschild black hole whereas the static 
exterior makes a positive contribution. Our numerical plots bear this out. The change in free 
energy in the interior region, the dip in fig. [H is negative and this occurs in the nonstationary 
region 3ft whereas in the exterior, which is basically static except for an outgoing matter wave, 
the change in free energy in fig. [J] is positive. We now compare numerical values of F = —L 
during gravitational collapse to the free energy of a Schwarzschild black hole. 

Figure [7] shows the time evolution of the function F = —L, along with its gravity and matter 
contribution. At the late time t = 20, F has a numerical value of 0.351. This needs to be 
compared to E/2, the value of the free energy of a AD Schwarzschild black hole of mass E. 
The value of E/2 is obtained directly from the horizon radius rg = E/2, the location where the 
lapse function N crosses at late times. Note how the plots of tq and F to t merge together at 
late times in fig. The value of rn at i = 20 is 0.334. The value F = 0.351 obtained numerically 
from — L is therefore equal to the E/2 result (0.334) to within 5%. The 5% discrepancy is a 
conservative value. It is based on a comparison between the free energy F and E/2 at t=20 
where Etot was just starting to deviate from unity. From t=20 to t=22, the numerical curve 
of the metric function tp in most of the exterior continues to approach closer to the analytical 
curve tp = 1 + M/2r. This suggests that the evolution from t=20 to t=22 is still good even 
though Etot has deviated from unity. In Fig. 7, the free energy F and the E/2 curve continue to 
approach each other in the time interval from t=20 to t=22 and the discrepancy at t=22 is 2.7% 
(F=0.344 and E/2=0.335)H. This can be viewed as a more realistic value for the percentage 
discrepancy. It is clear from this analysis that the main reason for the 5% discrepancy is that 
F and E/2 have not evolved far enough in time to reach a plateau before Etot starts to deviate 
from unity. The goal in the future would then be to push the time further numerically. It 
should be possible to accomplish this by increasing the resolution and implementing adaptive 
mesh refinement |25j to reduce the computing time (we ran the code with the larger step size 
of Ax = 0.001, At = 0.001 and the same A = 1.5. At t=20, the discrepancy is 14% compared 
to 5% and the E to t is quite off so that the higher resolution of Ax = 0.0001, At = 0.0001 leads 
clearly to better results). It is worth noting that the discrepancy also depends on the initial 
state. Recall that the initial field configuration for the scalar field x depends on the scale 
parameter A ()27p . We ran the code for different values of A (albeit at at the lower resolution 
of Ax = 0.001, At = 0.001 to speed up the simulation). We observed lower discrepancies for 
lower A in a given range (30% at A = 1.55, 12% at A = 1.5 and 10% at A = 1.49). The value 

3 The matter contribution to the free energy from t=20 to t=22 is separately evolving towards its expected 
theoretical value of zero: it is 0.003 at t=20 and 0.0004 at t=22. 
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of A also affects the duration of the simulation, extending the time where Etot is close to unity 
from T=15 for A = 1.55 up to T=17.5 for A = 1.49. 
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Figure 1: Space profile of the metric field ip for different times. At late times (t = 20), the 
exterior region matches almost exactly the Schwarzschild form. 
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Figure 2: Space profile of the lapse function iV for different times. The event horizon (r = 0.334) 
is the radius where iV crosses zero at late times. After crossing zero at r = 0.334, N is negative 
in the interior and approaches zero from below. 
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Figure 3: Space profile of the matter field x f° r different times. The initial matter distribution 
is a shell. It first expands into the interior (t = 4) before collapsing to a thin shell near the 
event horizon at late times. Note the propagation of an outgoing matter wave in the exterior 
region. 
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Figure 4: Accumulation plot for the gravitational Lagrangian at different times. At late times, 
there is a negative contribution (the large dip) stemming from a thin slice just inside the event 
horizon where the metric field tp is nonstationary. 
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Figure 5: Accumulation plots for the Klein-Gordon Lagrangian for different times. It is close 
to zero at late times. 
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Figure 6: Time dependence of the metric field ip. ip peaks just inside the event horizon, is 
decreasing in magnitude in the rest of the interior and is zero outside. 
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Figure 7: Time profile of the free energy and the horizon radius. The total free energy Ff t 
and the horizon radius r$ = E/2 converge towards the same value at late times. 

3 Free energy during the collapse of a 5D massless scalar field 

Though 3 + 1 dimensions represents the most realistic scenario, it would be short sighted not 
to consider the higher five-dimensional (4+1) case. Besides providing another arena where the 
thermodynamics of gravitational collapse can be studied numerically, 5D gravitational collapse 
has had surprising connections to 4L> physics. For example, a holographic argument backed 
up with numerical results was offered for a connection between the Choptuik critical scaling 
exponent, 75D, in 5-dimensional black hole formation via scalar field collapse |26} I27| and the 
saturation exponent, jbfkl, of four-dimensional Yang-Mills theory in the Regge limit [28] , 
The latest value of 75/5 = 0.4131 ± 0.0001 obtained from numerical work on 5D gravitational 
collapse [29j was close to the value of jbfkl = 0.409552 obtained from numerical work on 
AD Yang-Mills |28| . The near matching of the two results provides support for the conjecture 
relating 51? gravitational physics to pomeron exchange in 4D Yang- Mills theory [28]. We have 
learned by now that gravitational investigations in five dimensions are not just a curiosity but 
have something real to say about our four dimensional world. 

In this section we track the free energy F = —L during the classical gravitational collapse 
of a massless scalar field in five dimensional isotropic coordinates. Whereas the previously 
studied 5D Yang-Mills instanton [9] acts like dust and collapses to a black hole for any initial 
static configuration, a static massless scalar field has pressure and two outcomes are possible: 



IS 



dispersion or gravitational collapse to a black hole [13]. The initial configuration of the scalar 
field is governed by a scale parameter A and the system collapses to a black hole for sufficiently 
small values of A. At late stages of the collapse process, we obtain a numerical value of the 
function F = —L that is within 5% of the free energy of E/3 obtained from 5D standard black 
hole thermodynamics. Although a massless scalar field and a Yang-Mills instanton constitute 
different kinds of matter, they collapse to nearly the same final thermal state. We see in a 
dynamical fashion that black hole thermodynamics is not dependent on the type of matter 
undergoing the collapse. 

A spherically symmetric time-dependent 5D metric in isotropic coordinates takes the form 



where dVL 2 = d6 2 +sm 2 8(d(J) 2 +siyi 2 (pdj 2 ) is the metric on S s , the unit 3-sphere, with < 8 < 7r, 
< (j> < vr and < 7 < 2ir. As before, N(r,t) is called the lapse function and ip(r,t) the 
conformal factor. In 5D isotropic coordinates, the exterior region of the Schwarzschild metric 
is given by [9] 



where ro is the location of the event horizon and is given by rl = 2 G 5 £/(3tt) (G 5 is Newton's 
constant in 5D and E is the ADM mass of the black hole). The static metric (|38p does not 
describe the interior of the 5D Schwarzschild black hole, which is nonstationary. During the 
collapse process, the metric (|37|) approaches the metric (|38p only in the exterior region outside 
the event horizon. In the interior, the lapse function N(r, t) and conformal factor ip(r, t) deviate 
from (|38p and are time-dependent in accordance with the nonstationary nature of the interior. 

3.1 Gravitational and matter Lagrangian in 5D 

The Einstein-Hilbert action in 5D is given by 



ds 2 = -N(r, tfdt 2 + ij>(r, t) 2 (dr 2 + r 2 dnj) 



(37) 




(38) 




(39) 



where the Ricci scalar evaluated for the isotropic metric ([3?]) is 

R 5 = — ^-^r(Ar^ 2 N -4rNipi; 2 -9ip'N 3 -3ri()"N 3 (40) 



- 2rN'ij'N 2 - 3N'N 2 ip + 6rtp 2 N^ - rN"N 2 ip), 



and \f—g§ = Nip 4 " r 3 sin 2 9 sin (f>. Integrating by parts the second derivatives N" , ip" , and ip we 
obtain 
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where 

L G = JL. J ( - 12r W + 6r 3 iWV' + 6r 3 iV^' 2 ) (42) 

is the gravitational Lagrangian and Is is a boundary term [21]. Variation with respect to 
the metric of the gravitational action Iq = I eh + Ib = f Lq dt reproduces the Einstein field 
equations. 

The Klein-Gordon action for a massless scalar field in 5D is 

Ikg = ~\j ^gig^d^xduX^x = I ^g7 £ KG d 5 x (43) 

where spherical symmetry implies x = x( r it) anci £kg = (~ 1/2) 9^ duX^vX is the Lagrangian 
density. The matter Lagrangian evaluated for the isotropic metric (|37p is given by, 

Lkg = J V=9~5 Ckg d A x = -vr 2 J N ^ r 3 - ^) dr . (44) 
The total Lagrangian Lj D t = + Lkg is then 

L " = Wk J ( ~ "^JV 2 ~ + 6r3jy '^ + fo-W) ^ (45) 

During the collapse process we track the function F(t) = —Lt t- 
3.2 5D equations of motion 

In 4D, the Einstein field equations yield the energy constraint ([6]), the momentum constraint 
(fl"2|) . the evolution equation (fTD"]) for the conformal factor -0, the evolution equation ([To]) for its 
"momentum conjugate" , the trace of the extrinsic curvature K and finally the ODE (|17p for 
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the lapse function N. The corresponding equations in 5D are (see derivation in Appendix A) 

energy constraint: — 3 — 5- = k 2 £ K 2 (46) 

ip d 8 

K' k 2 x 1 

momentum constraint: — = ~^~~jyX (47) 

if) Kip 

evolution equation for tp : — = — (48) 

K K 2 4 
evolution equation for K: — = + w ) + 3^ ~ w ^ 

4 j// / jj' | 2\ ^ AT' /j// | 1\ iw Y^t (49) 
ip^\ip r) Nip 2 \ip r) ■0 3 

<ra fa N: (50, 

where the arbitrary parameter w in (|49p arises from a term added for numerical stability 
purposes and £ is interpreted as the energy density of the matter field and is given by 

Note that in 5D the flat space Laplacian is given by 

VV = ^3 Wr 3 ) = V + ^" ■ (52) 

There are also the evolution equations for the matter field x an d its "momentum conjugate" 
p. These are 

= 1 

N V 

and 



3.3 Initial states 



(53) 



N V[x" + ^x' + 2^ + -x') (54) 



The initial states for the matter and metric fields are chosen to be static: K=p = or tp = x = 
at t = 0. The initial configuration of the massless scalar field x is chosen to be a shell with 
Maxwellian distribution: 

x (r,t = 0) = ^r 2 e^ (55) 
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where the scale parameter A is a positive constant. The function peaks at r = V2 A and the 
half-width of the shell is ~ 1.6A. The function is zero at the two extremities: x( r > t = 0) — > as 
r — > and as r — > oo. A large A implies a diluted energy density and a weaker self-gravitational 
force. As in 4D, black hole formation occurs only for a sufficiently small value of A. 

The initial state for the conformal factor ip is obtained by substituting the above initial profile 
for x into the energy constraint (|46p and solving it order by order in powers of £ = k 2 as we 
did for the 4D case (we refer the reader to section 2.3 for the procedure and we will simply 
state the result here). The initial state for tp has to be expanded to at least order £ 3 for the 
simulation to be sufficiently accurate. Its expression is 



C -r 2 r 

i/j(r, t = 0) = 1 + ztAtth e * 



+ 
+ 



24 r 2 A 6 

c 2 

4608 r 2 A 12 

C 3 



-2r 

e ^ 



241864704 r 2 A 18 



. r 6_ 2r 4 A 2_ 6r 2 A 4 + 8A 6( e ^_ 1 > ) 

2r 10 + 6r 8 A 2 + . . . - 3r 2 A 8 + A 10 (77e^" - 128e^ + 5l) 

2869) 



e a 



486 r 14 + 1782 r 12 A 2 + . . . + 6 r 2 A 12 ( 11664 



- 2 A (105832 e^ 2 " - 168399e^" + 78732 - 16165) 

The initial state for N is obtained numerically from its ODE (|50p. using the initial states for 
ip and x an d iterating backwards starting with the asymptotic value N(r = R,t) = 1. 



3.4 ADM mass in 5D 

To obtain an expression for the ADM mass in 5D, one uses the same definition (|24j) but now 
integrated over d?9, 

Madm = —r- km f (k — k )^d 3 9, (56) 

87T 5t^oo J St 

where St is the three-sphere boundary. Let £t be the spacelike hypersurface at constant time. 
For the metric ([57]) . the trace of the extrinsic curvature of 5t embedded in St is given by 
A; = V a r a = 3tp' /ip 2 + 3/(ipr) where r a = , i/ ;< 9a^ is the unit normal vector to the three-sphere 
boundary St, while the trace of the extrinsic curvature of St embedded in flat spacetime is 
kg = 3/(tfir). With ^fo<$Q = tp 3 r 3 sin 2 9 shicpdO dcpdrf, the ADM mass for the 5D isotropic 
metric ([37]) reduces to 

MaDM = ~ dr ^r=R = -^Q- d Mr=R ( 57 ) 
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where tp = 1 at r = R was used. Using the energy constraint equation (|46p . one can express 
the ADM mass in the integral form 



The ADM mass represents the total energy E tot . As in 4D, the black hole mass E is less than 
Et t because part of the total energy is carried away by an outgoing matter wave. In 5D, E 
is obtained via the relation rg = 2G5 E/(3ir) where the gravitational radius tq is the location 
where the lapse function N crosses zero at late times. 

3.5 5D results 

We work in the following geometrized units: G^, = 1/(6tt) and c = l. In 5D geometrized units, 
mass has dimension of length squared (this is why in 5D, the black hole mass E is proportional 
to the square of gravitational radius tq). Masses will be expressed in terms of the ADM mass. 
The ADM mass depends on the parameter A which sets the length scale. We use the same 
numerical technique as in 4D (a fourth-order ABM scheme). We perform the spatial coordinate 
change r — > x = ^f2, which allows for more grid points near r = 1 where the gradients are 
largest. The space and time intervals used are Ax = 1 x 10 -4 and At = 5 x 10 -5 respectively. 
The value of the scale parameter A has to be small enough for gravitational collapse to occur 
and our simulation is run using A = 1 which leads to black hole formation instead of dispersion. 
The value of A affects the time of collapse but not the black hole thermodynamics at late times. 

Black hole formation occurs when the lapse function N in fig. [8] first crosses zero. The data 
shows that this occurs for the first time around t = 2.6. At late stages of the collapse, N crosses 
zero at r = 0.49, which is associated with the event horizon. In the interior region (r < 0.49), 
N becomes negative and tends to (but never reaches) zero. The spacetime is well behaved 
throughout the simulation. Figure [9] shows the time evolution of the matter field x- Early in 
the simulation, the shell expands into a spherical ball, filling the interior, before collapsing to 
a black hole at around t = 2.6. At late times, the matter field has collapsed to a thin shell 
immediately inside and on the horizon, except for an outgoing matter wave located at around 
r = 10 at t = 8. In the region outside the event horizon, the conformal factor ijj (fig. [T0|) 
matches closely the analytical Schwarzschild form given by metric (I37D so that the collapse 
process has indeed led to a Schwarzschild black hole. 

The temporal dependence of the metric field ip is shown in fig. QTJ At late times, the metric 
field tjj is nonstationary in a small region just inside the horizon. This reflects the fact that the 
Schwarzschild spacetime is nonstationary in the interior region. It is this small nonstationary 
region which is responsible for the large dip near the event horizon (r = 0.49) at late times 
in the accumulation plot of the free energy for gravity (see fig. \V2h . In accordance with our 




(58) 
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Figure 8: Space profile of the lapse function N at different times. N crosses zero at the location 
of the apparent horizon. This occurs for the first time at t = 2.6. At late times, N crosses zero 
at r = 0.49 and this is identified as the radius of the event horizon. After ./V crosses zero at 
r = 0.49, it is negative (albeit very small) throughout the interior (r < 0.49) and never crosses 
zero inside. 

discussion in section 2.4.1, the nonstationary interior region makes a negative contribution to 
the free energy, whereas the basically static exterior makes a positive contribution. The matter 
contribution to the free energy (see fig. fl~3l) is essentially zero except for a small contribution 
from the outgoing matter wave. The free energy is almost entirely gravitational in origin and in 
fig- US the total free energy Ftot is equal to the gravitational contribution Fq at late times i.e. 
the red and green plots meet at late times. Black hole entropy is gravitational and originates 
from the nonstationary interior region (see [2] for a further discussion on this point). 

The function F = —L at late times in fig. Q3]is equal to F = 0.1764. This needs to be compared 
to the free energy E/3 of a 5D Schwarzschild black hole of mass E, where E is calculated 
via the relation E = 3ir Tq/(2 G5) = 9tt 2 Tq (this is then divided by E to t to obtain its value in 
ADM mass units). E is plotted in fig. [14] and at late times (t = 8) its value is 0.5044 so that 
E/3 = 0.1682. There is less than a 5% difference between the numerical value of F = —L at 
late times and the expected E/3 result. 
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Figure 9: Space profile of the matter field x a t different times. At late times, the matter has 
collapsed to a thin shell near the event horizon. Note the propagation of an outgoing matter 
wave in the exterior region. 
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Figure 10: Space profile of the metric field ip at different times. At late times (t = 8), the 
exterior region matches almost exactly the Schwarzschild form. Note that tp — > in the 
interior. 
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Figure 11: Time dependence of the metric field ip. It is static except for a thin slice just inside 
the event horizon. 




Figure 12: Accumulation plot for the gravitational Lagrangian at different times. There is a 
negative contribution (the large dip) coming from the region just inside the event horizon. 
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Figure 13: Accumulation plot for the Klein-Gordon Lagrangian at different times. It ap- 
proaches zero at late times except for the small contribution of the outgoing matter wave. 




Figure 14: Time profile of the free energy and the mass of the black hole. The total free energy 
approaches a value of a third of the mass E at late times. 
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4 Conclusion 



In this paper, we tracked numerically the negative of the Lagrangian, —L. during the spherical 
collapse of a massless scalar field to a black hole in isotropic coordinates. This was performed in 
4D (3+1 dimensions) and in 5D and as far as we know, this is the first time that the physically 
relevant 3+1 case has been studied. We showed that at late stages of the collapse, — L reached 
a numerical value close to the free energy E/2 for the 4D Schwarzschild black hole (to within 
5%) and close to the free energy E/3 for the 5D case (again, to within 5%). This confirms 
that previous 5D thermodynamic results from the collapse of a Yang-Mills instanton [9] were 
universal and not dependent on the particular type of matter. We saw that the free energy 
stems almost entirely from the gravitational part — Lq of the total Lagrangian: the matter 
contribution, —Lkg, was close to zero. The main reason for the 5% discrepancy is that the 
fields have not evolved far enough in time before the total energy begins to deviate from unity. 
It should be possible to reduce the discrepancy in the future by increasing the resolution and 
implementing adaptive mesh refinement |25] to reduce simulation time. 

Our 4D and 5D thermodynamic results obtained from the dynamics of classical gravitational 
collapse together with previous 5D results [9] provide numerical evidence for associating a 
classical quantity, the negative of the gravitational Lagrangian —Lq, to the (Helmholtz) free 
energy F = E — T S of a black hole. This opens up an interesting question. Can the product 
TS have a classical interpretation without the temperature and entropy having a classical 
interpretation? It is sometimes stated that because black holes do not radiate classically, 
they must have zero temperature classically. However, it seems difficult to reconcile such a 
statement with what we know of the free energy i.e. that classically TS ^ 0. 

For a Schwarzschild black hole (vacuum solution), we showed analytically that the static exte- 
rior region makes a positive contribution E to —Lq. It follows that the negative contribution 
—TS to the free energy stems from the nonstationary interior region. This feature of the 
Schwarzschild black hole is observed in our numerical simulation. In our accumulation plot 
of —Lq, a negative contribution (a large dip) can be seen just inside the event horizon where 
the metric field ip is nonstationary. This ties in with Bekenstein's original view of black hole 
entropy as a measure of an outside observer's ignorance of the internal configurations hidden 
behind the event horizon [11]. As discussed in section 2.4.1, the internal configurations have 
recently been identified as points in phase space [h a b,P ab ] [H]. If the region behind the event 
horizon is nonstationary, then there is a continuous set of classical microstates [h a },(t) , P ab (t)] 
in the interior parametrized by the time t. An outside observer does not have access to the 
parameter t and is clearly ignorant of which classical microstate the black hole interior is in. 
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A Derivation of the 5D equations of motion 

In this appendix we derive the 5D equations of motion. To remain general, we do not set here 
the mass m of the scalar field to zero. To recover the equations of motion (|46p -(|50p simply set 
m = 0. 

We begin with the Einstein field equations 



The general expression for the stress-energy tensor T a p of the KG field is given by @ . From 
the Lagrangian density of the 5D KG action, along with the metric tensor of the metric (13?p . 
we obtain 



By calculating the components of this tensor explicitly, one can observe that only seven of 
them, two off-diagonal and the five diagonal elements, are non-zero: 



Gap — RaP — 77-^5 da/3 ~ ^TtG^T^ — K T a p. 



(A.l) 




(A.2) 




(A.3) 



(A.4) 



T, 



rt 



= T tr = ) 
■0 2 r 2 



XX 



(A.5) 




(A.6) 




sin 2 9 T ee , 
sin 2 9 sin 2 (f> Tee ■ 



(A.7) 
(A.8) 
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The components of the Einstein tensor are 

_3 

rip 



G a = -^{2ri> 2 4> -Zijj'N 2 -rif;"N 2 ), (A.9) 
3 



G rr = (-rib' 2 N 3 - 2ip'N 3 ib + rip 2 Nib 2 - rNipib 3 

rip z N d 

+ r^Nip 3 -rN'ip'N 2 ip - N'N 2 ip 2 ) , (A.10) 

G rt = G tr = ^-(-iPtP'N + ip l Ntlj-ipN'tP), (A.11) 

G ee = -^r(3rip 2 Nip 2 - 4ip'N 3 ip - 2rip"N 3 ip + rip' 2 N 3 - SrNipip 3 
ip z N 6 

+ 3ripNip 3 - rN'ip'N 2 ip - 2N'ip 2 N 2 - rN"N 2 ip 2 ) , (A.12) 

G H = sin 2 9Gge, (A. 13) 

G 77 = sin 2 9 sin 2 cpGee- ( A - 14 ) 



The (p, and 7 components of these two tensors can be expressed as a function of the 9 component 
and the rt and tr components are equal, which leaves us with only four independent equations 
(the ones for the tt, rr, rt and 99 components). 



A.l The energy constraint equation 

The energy constraint equation is obtained from the equation Gtt = k 2 Tu, and is given here 
as 

^(2nft& - 3^'iV 2 - riP"N 2 ) = ^{^ + f 2 + -V). (A.15) 

This equation can be reduced to a more convenient form by using the flat space Laplacian (|52j) 
and the trace of the extrinsic curvature in 5D isotropic coordinates obtained from the equation 

K = g ab K ab = ^g ab d t g ab = 

-1 /2ipip ^ 2ipipr 2 ^ 2ipipr 2 sin 2 9 ^ 2ipipr 2 sin 2 sin 2 <p \ 
2N \ ip 2 ip 2 r 2 tp 2 r 2 sin 2 9 ip 2 r 2 sin 2 9 sin 2 (p / ' 

which yields an evolution equation for the conformal factor ip 
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Dividing the energy constraint equation (|A. 15[) by iV 2 and using the Laplacian ()52[) and ()A.17|) , 
we obtain 

The term in parenthesis is, as in 4D, interpreted as an energy density of the matter field x 

e = \(w* + $ +v )- (A ' 19) 

where V = m 2 x 2 ■ 

The energy constraint equation in its final form is 

-^=JS-\ K \ (A.20) 

A. 2 The momentum constraint equation 

The momentum constraint equation is obtained from the Einstein equation G r t = 

3 (-ipip'N + ij/Ni/t-i)N'^) = k 2 xx' ■ (A.21) 



ip 2 N 

From the spatial derivative of the evolution equation for ip (|A.17|) we obtain 

K' _ / ip' N'tp tptp' 

T ~ ~ viv^ ~ N 2 ip ~ mp 2 



(A.22) 



By substituting the above equation in (|A.2ip . we obtain a more convenient form for the 
momentum constraint 

K 1 n 2 x i , . x 

T = (A - 23) 

A. 3 The evolution equation for K 

The evolution equation for K can be obtained by rearranging the Einstein equation G rr = 

~ 3 „ (rip 2 Nip 2 - 2ip'N 3 ip - rip' 2 N 3 - rNipip 3 + r^pNip 3 
rip z N 6 V 
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The time derivatives terms in the left-hand side can be expressed in terms of K and K from 
dQ7j) : 

tp 2 Ntpib 4>ip 3ip 2 (K K 2s 



-^ + ^- 3 ^ = i -(---j. (A.24) 

Substituting the above equation and the definition of £ (|A.19p in (|A.24p . one obtains an 
evolution equation for the trace of the extrinsic curvature K: 



— - — + -k 2 £ - -K 2 V - 4^ (— + -) - 4—2 {— + - 

N 2 3 3 ib 3 \ib r) Nib 2 \ib r 

From the prescription in [T7j, we add a term proportional to the energy constraint (|A.20|) to 
(IA.25P in order to get rid of numerical instabilities in the simulation: 

wC e = --w(k 2 £ 1 -r-), (A.26) 

3 3 V 8 V 3 

where w is a positive constant that has to be unity or very close to it. The final form of the 
evolution equation for K is then 

§ = ^{i + w ) + -^£(i- w )-\k 2 v 

A ib' fib' 2\ A N' 1\ A V 2 ip /A WN 

As = from (|A.20|) . the additional term will not alter the mathematical content of the 
equation, so ()A.25|) and (jA.27j) are equivalent in that respect. 



A. 4 An ODE for the lapse function 

An ODE for the lapse function N can be obtained from the last of Einstein's field equations, 
Gee = K 2 Tee, i.e. 

r -(rip' 2 N 3 - Aip'N 3 ip - 2rib"N 3 ib 



ib 2 N 3 

+ 3rip 2 Nib 2 - SrNipib 3 - rN'ib'N 2 ip 

+ 3ribNip 3 — 2N' 2 ip 2 N 2 — rN" N 2 ip 2 ) = - + m 2 x 2 ) . (A.28) 

2 V ib z N z J 

The combination of (|A.24p and (|A.25|) yields 

il 2 Nihil) ibib 

- 3 h +3 ^- 3 w= (A - 29) 

9 , ~ q \ fib' 2 2ib' N' ib' IN' 
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Substituting the above equation into (|A.28|) . one can obtain a convenient form for the ODE 
for N 

*d r (^-) + -d r (—) = -k 2 ^ 



M^.)+-A( — (A.30) 



A. 5 Equations of motion for the matter field 

The evolution equations related to the matter field x an d its "conjugate momentum" p are now 
the only equations left to derive. One way of obtaining them is from the Lagrange equation of 
motion 

9L t ot / dL to t 



Substituting the total Lagrangian ()46p and separating the spatial and time derivatives, this 
equation becomes 

4^- + " ^ = ^ [X" + W X' + 2^X' + -X' - ^x) ■ (A.32) 

The left-hand side of this equation can be expressed as 

Ax^ . X^ l a( ,4,X\ 1 a oo\ 

+ J^- — = N dt ^N)- (A - 33) 

A new function p is defined here as 

P = ^- (A.34) 

This function can be roughly associated to the momentum conjugate to x- (|A.32|) is then 
expressed by two single time derivative equations: 



X_ = V 
N ib 

and 



(A.35) 



P _„/.2/0' , N 'j, r^'j , 3 / ,2 2„ 



^ = r [x" + jj-x' + 2jx' + -X' - rm'x) • (A.36) 

Furthermore, the evolution equation for x (|A.35|) can be substituted in the energy density 
(|A.19|) and in the momentum constraint (IA,23j) to obtain 

1 (t- + ^ 

2 V</> 8 V> 2 

and 

4 3~0 4 



£ = ^( i Ts+^ + ^ (A.37) 



ft'' K 2 

^t = ^px'. (A.38) 
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